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Abstract 

Iterative solutions to fourth-order gravity describing static and electrically charged black holes are 
constructed. Obtained solutions are parametrized by two integration constants which are related 
to the electric charge and the exact location of the event horizon. Special emphasis is put on the 
extremal black holes. It is explicitly demonstrated that in the extremal limit, the exact location of 
the (degenerate) event horizon is given by r_|_ = |e|. Similarly to the classical Reissner-Nordstrom 
solution, the near-horizon geometry of the charged black holes in quadratic gravity, when expanded 
into the whole manifold, is simply that of Bertotti and Robinson. Similar considerations have been 
carried out for the boundary conditions of second type which employ the electric charge and the 
mass of the system as seen by a distant observer. The relations between results obtained within 
the framework of each method are briefly discussed. 
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I. INTRODUCTION 



According to our present understanding applicability of the conventional Einstein-Hilbert 
equations is limited to curvatures significantly less than the Planck scale and should be 
considered as the first approximation to a more fundamental theory. And although it is 
not clear how to construct this more fundamental theory, it seems reasonable to address 
the question of its possible low-energy impact. In the quest for imprints of the quantum 
gravity effects in the classical regime especially helpful observation is that regardless of the 
formulation of the fundamental theory its low energy effective action should consist of the 
classical gravity supplemented by covariant higher curvature terms and higher derivative 
terms involving other physical fields. The gravitational part of the total action can be 
written therefore as 

Ig = Io + h + I2 + (1) 

where /q is the Einstein-Hilbert action and (for i > 1) denotes combination of operators 
of dimension 2z + 2 with numerical coefficients depending on a type of the theory. 

Among the various modifications of general relativity proposed so far the prominent role 
is played by Lagrangians describing the quadratic gravity. The interest in the theories of this 
very type is motivated by the fact that such theories appear, as expected, in low energy limit 
of string theory and in the attempts to construct renormalizable theory of gravity coupled 
to matter lLy,yl. It should be also remembered that theories of this kind are almost as old 
a. ge„e.a. .e.atWitv ..e>f a.d a, *e .en... .e.ativ.y have a noUe parentage 3000. 

Higher derivative terms in the equations of the gravitational field have very important 
consequences. One of them consists in the obvious observation that such equations are pre- 
sumably hard to solve even in the cases when Einstein's equations admit solutions expressible 
in terms of known special functions. Further, it should be emphasized that although famify 
of solutions is richer one should be cautious in selecting physically meaningful solutions jsf. 
Naive acceptance of solutions may lead to interesting behavior of the system but a wrong 
physics. It is natural therefore that in order to gain insight into the nature of the nontrivial 
problem one has to refer to approximate methods. Such methods are expected to yield 
reasonable solutions to the problem simultaneously selecting the physical ones ^, ^, . 

For almost two decades perturbative methods have been extensively used in black hole 
physics in the context of the back reaction of the quantized fields (both massive and massless) 
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on the metric 
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13|,|l4,|li|l6|,|l7|,|l8|, 3- The (one-loop) approximation to the stress- 
energy tensor of the massive quantized fields in a large mass limit for example, could be 
regarded as higher curvature term constructed from the curvature tensor, its contractions 



and covariant derivatives 
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24l |. It should be noted, however, that the back 
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29|, and, 



reaction analyses are, inevitably, limited to the linearized equations. Similar considerations 
in the context of the string inspired action with quadratic and quartic terms contructed 
form the Riemann tensor have been carried out in Refs. I25I . |3| • 

The general perturbative method has been proposed in a series of papers 
subsequently, successfully applied in a number of physically interesting the de Sitter 

universe, cosmic strings, charged black holes, and gravitational waves js^. The method is 
simple: assuming that the quadratic terms could be considered as small contribution to the 
effective stress-energy tensor, which is justified for small curvatures of spacetime, one can 
iteratively solve resulting equations order by order starting with the classical Einstein field 
equations. Among various applications the most interesting from our point of view is the 
perturbative solution describing static, spherically-symmetric and electrically charged black 
hole constructed by Campanelli, Lousto and Audertsch (CLA) in Ref. 3]. The calculations 
have been carried up to third order with a special emphasis put on the thermodynamical 
issues. However, their analysis is somewhat obscured by the choice of the boundary condi- 
tions, and, moreover, propagation of errors in the metric tensor caused errors in some of the 
results. 

In this paper we shall return to this problem and express resultant solutions in terms 
of the electric charge and the exact location of the event horizon, r+. It seems that such a 
choice is more natural and reveals simplicity of the extremal configuration. It is also helpful 
in the analyses of the near-horizon geometry of the extremal black hole. Our method is 
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31j with the different choice of the boundary 



similar to that of York •11 and Lu and Wise 
conditions. 

The paper is organized as follows. In Sec. |n] we introduce basic equations and briefly 
sketch employed method. We choose the line element in the form propounded by Visser , 
which has proved to be a very useful representation. Since the functional dependence on the 
metric tensor of all terms appearing in equations is known, we start with the exact location 
of the event horizon from the very beginning. It means that the higher order terms do not 
contribute to r_|_. To express our results in a more familiar form we also introduce the horizon 
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defined mass. Discussion of tlie temperature and entropy of the nonextremal configuration 
is presented in Sec. IIIII A careful examination of the extremal configuration and the role 
played by the Bertotti-Robinson geometry is given in Sec. IIVI Moreover, it is shown that 
the extremal black holes are characterized by r+ = |e|. In Sec. |V]we discuss the problem 
from the point of view of a distant observer and explicitly demonstrate relations between 
both choices of the boundary conditions. Sec. IVII contains final remarks. 



II. EQUATIONS 



The coupled system of the electrodynamics and the quadratic gravity with the cosmolog- 
ical term set to zero is described by the action 



IGvr 



where 



and 



g'/^Fd'x. 



(2) 



(3) 



(4) 



Here F = FabF"-'^, Fab = ^b,a ~ ^a,b and all symbols have their usual meaning. The 
Kretschmann scalar has been relegated from the action by means of the Gauss-Bonnet 
invariant. 

Of numerical parameters a and (5 we assume that they are small and of comparable 
order otherwise they would lead to observational effects within our solar system. Their 
ultimate values should be determined from observations of light deflection, binary pulsars, 
and cosmological data 0, 33, Is^. Following [2^ we shall restrict ourselves to spacetimes 
of small curvatures, for which conditions 



|a;i?| << 1, 



\(iRab\ « 1, 



(5) 



hold. Additional constraints could be obtained from a non-tachyon conditions, which are 



closely connected with stability of solutions 



possess a real mass 
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. Demanding the linearized equations to 



37|, one obtains 



3a - /5 > 0, /3 < 0. 



(6) 



To simplify calculations, especially to keep control of order of terms in complicated series 
expansions, we shall introduce another (dimensionless) parameter e substituting a ^ ea 
and (3 ^ e p. We shall put e = 1 in the final stage of calculations. As the coefficient a 
does not appear in the final formulas, introduction of the additional parameter might appear 
as an unnecessary complication in the present context. However, it is really helpful when 
dealing with more general Lagrangians, as, for example, these of nonlinear electrodynamics. 

Differentiating functionally S with respect to the metric tensor one obtains the system 



where 



Gl-a^'^Hl- P^'^Hl = S7:Tl (7) 



j g^'^Rd^x = -G"^^ (8) 



2R' - 2RR"'^ + - c/"^ {R^ - AUR) (9) 



and 



5 ' 5goh J 



^; oh _ g^afe _ 2R^^R''^'^ + 1^°^ {RcdR"'^ - UR) . (10) 



The Faraday tensor Fab satisfies 



= (11) 



and 



*F"';a = 0. (12) 

The stress-energy tensor, T"-^, defined as 

T — — JT^T Sem (13) 

is therefore given by 

Let us consider the spherically-symmetric and static geometry. As is well-known the 
spacetime metric can be cast in the form 

ds^ = -/ (r) e^^^'-^de + ^ + {de^ + sin^ ^^0^) (15) 



/M = i-^. (16) 

The metric has horizons at values of the radial coordinate satisfying 

m(r+) = ^. (17) 

In what follows we shall restrict our analyzes to the outermost horizon, i. e., the one for 
which 

T 

m{r) < - (18) 

for r > r+. 

Spherical symmetry places restrictions on the form of the tensor Fab'- the only nonvan- 
ishing components of the Faraday tensor are connected with the static radial electric and 
magnetic fields. Simple integration of the Maxwell equations gives 



^01 



.2 



(19) 



and 

F23 = 02 sin ^, (20) 

where ai and 02 are integration constants interpreted as the electric and the magnetic charge, 
respectively, and the former will be henceforth denoted by e. In what follows we shall confine 
ourselves to the solutions with an electric charge only. 

The stress-energy tensor for the line element ()15|) and the Faraday tensor (|T^ is simply 

T/ = t; = -T,' = -r; = (21) 

and in this form is independent of the metric potentials. 

Now we are going to construct approximate solution describing static and spherically- 
symmetric charged black hole. The system of differential equations for m(r) and '0(r) is 
to be supplemented with the appropriate, physically motivated boundary conditions. Our 
preferred choice, requiring knowledge of the exact location of the event horizon, r_|_, is 

m(r+) = ^, (22) 

which, as we shall see, is related to a horizon defined mass of the black hole. With such a 
choice we express the solution in terms of the exact location of the event horizon and the 



electric charge. On the other hand it seems natural to express the results in term of the 
mass of the system as seen by a distant observer: 



lim m(r) = M^. (23) 

r— >oo 

For the function '0(r) we shall adopt a natural condition 

lim i){r) = 0. (24) 

r^oo 

Let us select the boundary conditions of the first type. Of functions m(r) and if){r) we 
assume that they could be expanded as 



m(r) = Mo(r) + ^£'=Mfc(r) + C(e"+^) (25) 

k=l 

and 

n 

V^(r) = ^£^^,(r)+ 0(5"+i) (26) 
fe=i 

It should be noted that ipo{r) = 0. 

Since we have assumed the expansion of m(r) in the form given by eq. (|25p. the condition 
()22|) could be rewritten in the form: 

/-^ if^ = 0, 

M,(r+) = { ' (27) 
if i > 1 



and such a choice is a typical mathematical procedure Now, inserting the line element 
flT^ into Eqs. ((Tj), making use of the expansions and and finally collecting the 
terms with the like powers of the parameter e, one obtains the system of the differential 
equations for Mq, Mk, and ijjk {k > 1) of ascending complexity. The zeroth-order equations 
reduce to that of Einstein-Maxwell system, and, after simple manipulations, lead to the 
following integral 

1 ^2 



Mo = - I r'Fdr + Ci, (28) 



which, when combined with the condition ()22|1 . yields 
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The first-order equation constructed from tlie radial component of Eq. (|7j) reads 
_ J2M'^ SMqMI, , 2M^^ 2K , 5MoK KK 

— pi •\ ^ 2 ~'~ 2 

I lyiLi rytj ry^ ly ry 

+ Mt + m(^) - ^^^^ - Mf ) + r m(^) - 2 Mo M^^^ j . (30) 

On the other hand, the first order equation constructed from the time component of Eq. ((7j) 
when combined with (pUj). could be easily integrated to yield 



Mr) =/?(Mf)-lM^) (31) 



where C2 is an integration constant. 

Inserting Eq. into (jHUj) and integrating the thus obtained equation one has 

^/2e2 3eV+ 3 Ge^ 3eM 



^ r3 2r4 2r+r4 Sr^ 2r^ lOr^^ 
Further, substituting the zeroth-order solution to Eq. (j3T| and making use of the boundary 
condition (f^^ gives 

Mr) = f3^,- (33) 

Starting from the second-order the differential equations become more and more complicated 
and we shall display their solutions only. After some algebra one has 

^2 /^351 eV+ 36 1156 76 eV 76 704 , eV^ 

^ ~ r+ r8 ^ 2^^ ~ 10 ^ 12 ^ ~ 28 ^ 

and 

^2 = (32^ + 32^ - 24^ - 4l4V (35) 
for the second order, whereas the third order results read 
_ , , f eV+ 12392664 331624 2229 



55 220 j^Jj.^ j.b j,-^ 5 ^6^^ 

..^o 130732 eV, 130732 115176 eV^ 51347 e^r^ 

5508 r -f -7^ h 



115176 51347 eVl 
+ ^ ^ -r - 2744 ^— + 80 + 32 ^- + 6816 



^ ]_ 1 ^ "^"1" 'y^^ ryi ly^lO ry'J ry^ lyT ryQ 

351 1053 9 \ 
- 48 ^ H ^ —5- H ^ (36) 

3 ^7 4 y,3 ^8 20 rir^ 4 r'^ri 28 rlr^ / 



and 



, , , ,'96 , eV, 74560 eV+ 74560 

V'sW = ——^- 3264— ±- 3264- 32 



e^rl 48384 69572 e^X 

+ 2352-^ + ^ — + 2352-^ + -— — + 1080-^ . (37) 



It is possible to express this result in a more familiar form introducing the horizon defined 
mass M, i. e. to represent solution in terms of (e, M) rather than (e, r^). It could be easily 
done employing equality 

M ^ !± + ^ (38) 



and repeating calculations with 



Mo(r) = M - (39) 



Now the exact location of the event horizon is related to the horizon defined mass by a 
classical formula 

r+ = M + (M^ - e^f^. (40) 

It should be noted, however, that although the zeroth-order equation gives the exact location 
of the event horizon, the same is not true for its second root, 

^2 



e 



(41) 



Indeed, it could be shown that the inner horizon, r_, is given by 

(3 /3e^ 2e2 2 2r+ ?>rl\ ( 214 

r_ = r, + 3 ± + _i: +/52 _ 

5 \ r^J. r\ r+ e'^ / \ br\ 525 

38 206 131 58 r+ 374 542 191 

J \ \ ± J ± J ± + + Oid ) 

525 rl b2br\ 105eV+ 105 525 525 150 e^o ) ^' 

(42) 

and, moreover, simple calculation shows that for r+ = |e| both horizons coincide to required 
order. Such a behaviour strongly suggests that this very relation describes degenerate hori- 
zon of the extreme black hole. Since the extreme black holes deserve more accurate treatment 
we shall postpone further analysis of this and related problems to Sec. |TVl 



9 



III. TEMPERATURE AND ENTROPY 



One of the most important characteristics of the black hole is its Hawking temperature, 
Tfj. To investigate how Th is modified by quadratic terms, we employ the standard procedure 
of the Wick rotation. The thus obtained complexified line element has no conical singularity 
as r r+ provided the time coordinate is periodic with a period P, which is to be identified 
with the reciprocal of T^- Elementary considerations carried out for the general static and 
spherically symmetric spacetime describing a black hole lead to the formula 



For the line element (fT3jl and the boundary condition ()22|) the above result assumes simple 
form 



Making use of Eqs. and (jS2HSZI), and collecting the terms with the like powers of j3, 
after massive simplifications, one obtains 



We suspect that the common factor in the above expression appears in all higher-order 
terms and for r+ = |e| the black hole temperature approaches zero. We shall return to this 
problem later in the text. 

We have not, as yet, imposed any constraints on the parameter (3, but now we are going 
to examine the consequences of the non-tachyon condition, which, in the case at hand, is 
simply /? < 0. For the nonextremal black hole the terms proportional to (] and in the 
right hand side of the equation are strictly negative, and, therefore, one concludes that 
the temperature (to this order) is lower as compared with the Reissner-Nordstrom black 
hole described by the same values of r+ and e. On the other hand the higher-order terms 
proportional to (3'', (for A; > 3) change sign, however, for reasonable values of the coupling 
constant their contribution to a total temperature is negligible. 

The entropy of the black hole in quadratic gravity may be calculated using various meth- 
ods. Here we shall employ the Wald's Noether charge technique js^, which, as has been 




(43) 




(44) 




(45) 
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shown in Refs. 



40 



4l|, may be safely applied for Lagrangians of the type and leads to 



remarkably simple and elegant general result 

S = \ f d^xVh [1 + 2aR + (3 [R - If^Rab)] , (46) 

where hab is induced metric on S and the surface integral is taken across arbitrary section of 
the event horizon. In our representation, the entropy, after some algebra may be compactly 
written as 

^ ^ ^ _ 8^^, _ 1024^ . 

A A A^ ^ ' 

+ ^^'^^^''^f^' (A _ 4716^)2 (267re^ - 5A) + 0(/?^), (47) 

where A = Anr'^ is the surface of the event horizon. We attribute discrepancies between 
our result and the entropy computed by CLA to errors in their metric tensor. 



IV. THE EXTREMAL CONFIGURATION AND ITS NEAR HORIZON GEOME- 
TRY 

In this section we shall investigate the important issue of the extremal black holes. On 

general grounds one expects that the extremal configuration is the one in which (at least) 

two horizons merge. It means that in the simplest case the first two terms of the expansion 

df A d'^f 

f{r) = /(r+) + -f (r - r+) + (r - r^)' + ... (48) 

vanish. It is evident that if first n terms in the above expansion are absent one has a fold 
merging of n horizons. In a view of the further applications we also expect regularity of the 
function ilj{r) as the degenerate event horizon is approached: 

|V'(r+)|<oo, 1^ |<oo. (49) 

CLT |r=r_|_ 

In order to determine the location of the event horizon, i. e. to relate integration constants 
r_|_ and e one has to consider consequences of vanishing of the black hole temperature (surface 
gravity). Eq. strongly indicates that r+ = |e|. To demonstrate validity of this relation 
let us consider the (rr) component of ((Tj). It could be shown that if /'(r+) = 0, the conditions 
(j49|) are satisfied, and, additionally, 

|/(3)(r+)| < oo, |^(^n^+)l < oo, (50) 
11 



then the second derivative of / computed at r = r+ satisfies the constraint equation 



A \ dr"^ |r+ 



(51) 



The metric of the electrically charged black hole is to be independent of the coupling constant 
a 13, 0|- Repeating calculations for the (rr) component of the tensor ^^''H^^ at r = r+ one 
has 

1 rd'f ^ 



(52) 



2 \dr'^ \rj, 

Absence of terms proportional to {k = 0,1,2...) requires ^^^H^ to vanish identically, 
and, since the both conditions are to be satisfied simultaneously one obtains expected exact 
result: 

r+ = |e|. (53) 
Equivalently, expressing the location of in terms of the horizon defined mass yields 



M. 



(54) 



Obtained relations are identical with those which characterize geometry of the extremal 
Reissner-Nordstrom. As the particular form of the traceless stress-energy tensor played the 
decisive role in our derivation one should not expect simple generalizations. In a more general 
case, as for example that of the nonlinear electrodynamics, there is an explicit dependence 
on the coupling constant a, and the stress-energy tensor lacks simplicity of it Maxwellian 
analog. Indeed, repeating calculations for a line element of the type (fT^ and an arbitrary 
stress-energy tensor, one has 



and 



iP + 2a) 



- (/? + 2a) 



1 \ 1 



1 



(55) 



(56) 



It should be noted that in order to obtain Eq. ()5fij) one has to assume |/*-'^''(?"+)| < oo 
and |'0^^''(r_|_)| < oo. The spherical symmetry and the field equations at the event horizon 
respectively give 

T,' = Ti, T/ = r;. (57) 
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Combining the above equations one obtains 

- - = 87rr:(r+), (58) 

where is the trace of the stress-energy tensor. 

As is well-known the closest vicinity of the event horizon of the extremal Reissner- 
Nordstrom black hole, after the coordinate transformation 

r = r^[l + '-±y (59) 



could be approximated by the Bertotti-Robinson line element 



ds" = 4 i-de + dy^ + y^dn'') . (60) 

y 

Employing new static coordinates the Bertotti-Robinson line element can be rewritten in 
the form 

ds^ = rl (- sinh^ xdr^ + dx^ + d9^ + sin^ Odcf)'^) (61) 



or 

-1 



ds^ = -(^^-1^ dT^ + dx' + rl {de' + sin' ed<p') . 



(62) 



The latter form is particularly useful in demonstrating that topology of the Bertotti- 
Robinson solution is AdS2 x S^, i.e., it is a simple topological product of a (1 + 1) -dimensional 
anti-deSitter spacetime and two-sphere of radius r+. 

One expects that this spacetime also plays some role for the extremal black holes in the 
quadratic gravity. Indeed, first observe that because of simplicity of the Bertotti-Robinson 
line element the tensors ^'^^Hab and ^'^^Hab as well as the curvature scalar vanish. The only 
contribution to the left hand side of ((Tj) comes, therefore, from the Ricci tensor, which for 
the metric ()60|) is given by 

= K = -Re = -R^ = — 2- (6^) 

As the stress-energy tensor of the electromagnetic field is simply 



n = T; = -T,' = = (64) 

one concludes that the line element ()60|) is an exact solution if r+ = |e| . Moreover, inspection 
of ^ suggests that, after redefinition of the time coordinate, the dominant contribution to 

13 



the metric potentials gu and grr near the event horizon comes from the function Mo(r). 
It follows then that locally it resembles the geometry of the extremal Reissner-Nordstrom 
black hole. 

To investigate the role played by the Bertotti-Robinson solution in more details let us 
return to the expansion (|48|) and ask when the near horizon geometry of the extremal black 
hole is that described by the line element ()6U|). In the vicinity of the event horizon the line 
element may be written as 

ds^ = -e^^ir+)F{r - r+fdf + -dr^ + rUn^ (65) 

F[r — r^y 

where 

It could be easily demonstrated that expressing the line element in terms of a new coordinate 
y defined by means of the relation 



r = r I 



one obtains the line element ()60|) provided 

dr'^ \r=r+ r\' 

Inspection of ()52j) shows that it is precisely the relation which are satisfied by extremal 
black hole. It should be noted that by ()58|) our demonstration requires the trace of the 
stress-energy tensor to vanish at r = only. 

On the other hand, if does not obey Eq. (jUH)) . the topology of the solution (jU^ 

is still a simple product of AdS2 x S^, but with different modulus of curvatures. Indeed, it 
could be easily shown that 

R = i^AdS2 + Ks2, (69) 

where 

i^Ads, = -2F and = ^. (70) 

r_,_ 

Upon substitution 



one obtains 



ds"" = {-dt^ + dy^) + rl {de^ + sin^ 9d(/)^) . (72) 
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Finally observe that vanishing of the curvature scalar as r — > r_|_ yields ()68p . and the line 
element f|72|l reduces to (pH). 



V. A DISTANT OBSERVER POINT OF VIEW 



In this section we shall briefly examine consequences of the second choice of the boundary 
conditions as given by Eqs. ()23|) and ()24|) . This problem has been considered earlier by CLA. 
Unfortunately, their metric tensor and hence other characteristics of the charged black holes 
in quadratic gravity contain errors in the terms proportional to P^. 

Before proceeding further let us return to the function m(r) constructed with the aid of 
the condition ()22j). Its limit as r — oo is simply 



~ ~ 2 ^ 2r+ ^ ^[lOrl 2rl) ^ ^ \l2rl 28rl 
/4e^ _ 652e^ 2587e^ _ 2229e«\ 
^ ^ rl 55rl + 220rii 572r? ) + ^^^^ ^' i^^^ 



and is interpreted as a total mass of a system as seen by a distant observer and expressed in 
terms of the exact location of the event horizon and electric charge. On the other hand one 
can represent solution to the system ((7j) in terms of the total mass as seen from large distances 
and the electric charge from the very beginning. In this case the boundary conditions for 
the expansion ()25|) are to be rewirtten in the form: 

Moo if « = 0, 

lim Mi{r) = { (74) 
if i > 1, 



r— s-oo 



while the condition (j24j) remains, of course, intact. Repeating calculations order by order 
with the new boundary conditions one obtains 



m(r) = M^-- + P 2- -3^^+ ^ 
, .2 A ,0 596 704 e'Mj 351 e'^M^ 



+/33 ( 27264 ^^7^ - 11016 + 1440 % + 



4330344 179144 e^M^ 



\ J.7 3gg g ^10 

331624 460704 Moo ^ 51347 Moo e^Moo^ „ ,^,^4, 

+ ^ - ^ ^ - 21952 ^ + 7536 - j + 

(75) 
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and 



^(0 = P-i - \ —p. — + 



ly'l lyS 



/24864e^ 9408e2M2 695726^ 



5j,io ^10 ISr-*^^ 

1080e^ _ 6528e^Moo _ U9120e^M^ 



(76) 



Eqs. (f7S|) and (fTUj) are sufficient to determine ^fft and (^^r to The thus obtained metric 

tensor differs from that obtained by CLA. To demonstrate that our calculations lead to 
correct results let us make use of the consistency check. Inserting Eq. (fTSj) into (f75j) and 
fl7H|) one obtains, as expected, Eqs. and (jH^HTfj) . Since the calculations have been 
carried out independently we conclude that our results are correct. 

To determine location of the event horizon one can either iteratively solve the equation 

9u{r+) = (77) 

or revert relation ()7Hj) . Both methods give the same result, which reads 

^ (5r2 - 3e2) (2925r6 - 6515ry + 5095r2e^ - 13376^) 



5r|] (r2 - e2) ^ iQSOr^ (r^ - e^f 

/ 17268913 e^^ 61234643 e^^ 137993 



e2)^ V 75075 750750 r-j 770 tq 



3265043 26686967 436497 e^^ , 3064 2 , .^,^4 

+ ^ ^ + TT - 8 r;^ eVo + C 

10010 75075 35750 r^^ ° 55 ' 



(78) 



where 



ro = Moo + (M^-e^)'/'. (79) 

Repeating the steps of Sec. 11111 necessary to compute the Hawking temperature, and 
expressing the final result in terms of e and tq, one has 

(2r2 - e2) 



47rrQ ' 5rQ (rg 

7,^11 (^2 _ g2)3 ^25 21 ^-^' 75^-'^ 525 ■^■'^ 7 



(80) 
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Although the extremal configuration can be studied in (e, Moo) representation it is not the 
best choice. Indeed, simple calculations give for the location of the event horizon 

1 17 317 

whereas relation between the total mass as seen by a distant observer and the electric charge 
has the form: 



This could be contrasted with the simple relation between r+, and |e| given by Eqs. (j^H)) 
and ()54|). Finally observe, that since the entropy as given by Eq. ()47p is expressed in terms 
of the surface of the event horizon and the electric charge, we expect that it should be 
described by the same formula for both choices of the conditions (j^^ and (j^Hj) . 



VI. CONCLUSION AND SUMMARY 



In this paper we have constructed iterative solutions describing spherically-symmetric 
and static black holes to the equations of the fourth-order gravity with the source term 
given by the stress-energy tensor of the electromagnetic field . Obtained line elements 
are parametrized by two integration constants which are related to the electric charge and 
the exact location of the event horizon. The thus computed metric potentials enabled 
construction of the basic characteristics of the black hole: its Hawking temperature and 
entropy. 

Special emphasis has been put on the extremal black holes. Specifically, it has been ex- 
plicitly demonstrated that in the extremal limit, the exact location of the (degenerate) event 
horizon is given by r^. = |e| . It was shown that, similarly to the classical Reissner- Nordstrom 
solution, the near-horizon geometry of the charged black holes in quadratic gravity, when 
expanded into the whole manifold, is simply that of Bertotti and Robinson. As a byproduct 
of our investigations we obtained a simple equation that relate horizon value /" with the 
trace of the stress-energy tensor, which, for solutions of the equations of quadratic grav- 
ity describing black holes, may serve as an useful criterion for possessing the near-horizon 
geometry of the Bertotti-Robinson type. 

Similar considerations have been carried out for the boundary conditions of second type 
which employ the electric charge and the mass of the system as seen by a distant observer. 
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Moreover, it has been explicitly demonstrated how to relate appropriate results obtained 
within the framework of each method. 

Returning to the extremal black holes we observe that there are good reasons to believe 
that they are qualitatively different from the nonextremal ones. A proper distance, for 
example, between two points, one of which resides on the event horizon, is infinite. It could 
be easily seen from the integral 

J ^/g^dr ~ -^ln(r-r+), (83) 

as it diverges in the limit r — > r_|_. Moreover, the entropy remains nonzero as A — > Anr"^, 
and depends on the electric charge. This behaviour clearly violates the Nernst formulation 
of the third law of thermodynamics, which states that the entropy of a system must go 
to zero or a universal constant as T ^ 0. On the other hand, even if the entropy of the 
extremal black hole vanish 45] there are still problems simply because of its noncontinous 
nature. Indeed, zero is not the limit to which the entropy of the nearly extremal black holes 
tends. One can argue that this behaviour should not be treated as to worrisome simply 
because the Nernst formulation probably should not be considered as the fundamental law 
of thermodynamics. For a recent discussion of this issue see [4^]. On the other hand 
according to more radical opinions this failure indicates that the extremal black holes must 
not be treated as thermodynamic systems to which one can assign the notion of temperature 
and entropy |47|. Of course, we are unable here to judge which option should be treated 
seriously. All we can say now is that the entropy as given by (jTTj) (probably) could be safely 
used for nonextremal black holes. However, the particular case of the extreme black holes 
certainly deserves further study. 

Finally, let us observe that the methods of this paper could by easily generalized to other 
sources, such as, for example, nonlinear electrodynamics. This group of problems is under 
active investigations and the results will be published elsewhere. 
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